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SUMMARY 

One  of  the  important  variables  to  be  determined  when 
monitoring  unknown  radio  transmissions  Is  the  modulation 
type.  Automatic  recognition  procedures  based  on  time 
domain  parameters  additional  to  the  standard  envelope  and 
instantaneous  frequency  are  proposed.  Initial  simulations 
which  incorporate  Gaussian  noise  perturbation  and  centre 
frequency  and  bandwidth  mismatch  support  the  theoretical 
arguments  that  the  new  parameters  provide  better 
modulation  type  discrimination  in  noisy  conditions. 
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1 .  INTRODUCTION 

Vfhen  unknown  radio  signals  are  monitored,  either  for  military  purposes  or  to 
enforce  civilian  compliance,  one  of  the  important  variables  of  the  trans¬ 
mission  is  modulation  type(ref.l).  The  problem  of  replacing  operator 
intervention  by  machine  recognition  has  recently  received  attention.  Pattern 
recognition  procedures  based  on  time-averaged  behaviour  of  instantaneous 
envelope,  frequency  and  zero-phase  have  been  successfully  applied  to  strong 
signals(ref . 2  to  4). 

Existing  recognition  procedures  are  less  successful  at  low  levels  of 
signal-to-noise  ratio  (SNR) ,  in  part  because  these  methods  depend  on 
parameters  which,  in  the  presence  of  bandlimited  Gaussian  noise,  are  biased 
estimators  of  the  true  signal  parameters.  Moreover,  the  standard  approach  to 
modulation  recognition  is  inflexible  in  its  explicit  dealing  with  SNR. 
Typically,  the  signal  envelope,  frequency  and  zero  phase  are  sampled  at  the 
Nyquist  rate  for  a  total  duration  of  say,  0.1  to  1  s.  The  mean,  standard 
deviation  and,  possibly,  digitised  distribution  of  each  of  these 
characteristics  are  computed.  These  values  form  a  "feature  vector"  of  the 
signal,  which  is  compared  with  standard  feature  vectors  of  each  of  the 
modulation  types  being  considered  eg  FM,  AM,  SSB,  DSB,  FSK.  The  classifi¬ 
cation  of  the  signal  is  made  according  to  some  minimum  distance  criterion 
after  some  transformation  of  the  space  which  attempts  to  weight  components  of 
the  vector  according  to  their  discriminatory  value.  This  ignores  the 
intuitively  obvious  point  that  noisy  signals  are  more  similar  to  each  other, 
regardless  of  modulation  type,  than  they  are  to  strong  signals  of  the  same 
modulation  type. 

In  this  preliminary  report,  we  propose  the  use  of  a  modulation  recognition 
procedure  based  on  time-domain  signal  parameters  which  are  unbiased  estimators 
of  the  true  signal  parameters  in  the  presence  of  Gaussian  noise  with  symmetric 
spectral  density.  Time  averages  and  digitised  distributions  of  these  para¬ 
meters  would  be  used  to  form  feature  vectors,  as  outlined  above.  However,  we 
propose  assigning  to  each  modulation  type  multiple  classes,  corresponding  to 
various  signal  strengths. 

In  Section  2,  we  define  the  standard  parameters  and  indicate  theoretically  why 
features  based  on  them  are  unreliable  discriminators  of  modulation  type.  In 
Section  3,  we  introduce  some  unbiased  time  domain  parameters.  In  Section  4  we 
set  up  various  modulation  types,  and  in  Section  5  describe  the  signals  used  in 
the  limited  simulations  done  so  far.  Results  of  those  simulations,  presented 
in  Section  6,  support  the  theoretical  indications  that  the  new  approach 
enhances  recognition  of  modulation  types. 


2.  SIGNAL  CHARACTERISTICS  AND  SNR 
Let  the  received  signal  r(t)  be 


r(t)  -  A(t)cos6(t) )  =  A(t)cos(0' (t)t  +  6(t)).  (1) 


The  decomposition  of  equation  (1)  can  be  made  for  any  real-valued  function  r 
by  defining  the  envelope  A(t)  and  phase  0(t)  or  instantaneous  angular 

A 

frequency  9'(t)  and  zero  phase  6(t)  in  terms  of  r,  its  Hilbert  transform  r; 

A 

and  the  time  derivatives  r'  and  r',  viz: 
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A(t) 

-  (r(t)2  +  r(t)2)‘ 

0(0 

A 

-  arctan(r(t) /r(t) ) 

e'(t) 

-  (r(t)r'(t)  -  r (t) r ' (t) ) /A2 (t) 

6(0 

-  e ( t )  -  e* (t)t. 

(2) 

For  narrowband  signals  on  a  relatively  high  frequency  carrier,  the  Hilbert 
transform  can  be  obtained  instantaneously  by  analog  means  as  in  figure  1.  In 
this  case,  the  envelope  and  instantaneous  frequency  computed  using 
equation  (2)  coincide  with  these  signal  characteristics  as  obtained  directly 
by  analog  devices. 


Figure  1 .  Analog  computation  of  the  Hilbert  transform  pair 

Now  suppose  the  received  signal  r(t)  consists  of  a  modulated  signal  S(t)  with 
envelope  B(t)  and  phase  $(t),  plus  a  narrowband  Gaussian  noise  n(t)  of  power 
c2.  Then  it  is  well  known  that  the  expected  value,  EfA(t)),  of  the  envelope 
A(t)  at  any  instant  t  is 

E(A(t) )  =  o  M(-!,l,-B2(t)/2o2)  (3) 


where  M  is  the  confluent  hypergeometric  function.  (See  for  example 
reference  5,  page  105).  It  is  straightforward  to  show  that  E(A(t))  increases 
monotonically  with  a2,  having  asymptotic  expressions  B(t)  +  0(2a2/B(t)) 
when  a2  <<  B2(t)  and  /( to 2 / 2 )  (1  +  Bz(t)/4o2)  +  0(B‘*(t)/32a3)  when 

2k 

a2  »  B2 (t) .  Similarly,  for  k  >  1,  E(A  (t))  increases  monotonically  with  a2, 

2k  k 

from  B  (t)  at  a2  ”  0  to  the  asymptote  k ! ( 2o 2 )  •  Specifically, 

E (A2 ( t ) )  -  B2 ( t)  +  2a2 

and  (4) 


E(A**(t)) 


B“(t)  +  8o2B2(t)  +  Bau • 


■ST 


-  3  - 
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Rice's  classic  analysis  of  FM  modulation  showed  that  the  instantaneous  fre¬ 
quency  is  a  biased  estimator  of  the  signal  parameter  in  the  presence  of 
narrowband  Gaussian  noise  with  power  spectrum  symmetrically  distributed  about 
the  centre  frequency  fc(ref.6): 


E(6'  (t)-2irf  ) 
c 


(1 


-B2(t)/2a2 

e 


)(<*>' (t)-27Tfc). 


(5) 


To  get  a  feel  for  the  distortion  involved,  consider  a  signal  received  with  a 
SNR  of  0  dB  ie  B(t)  =  i'2c2.  Then  the  expected  value  of  0 *  ( t )  — 2rrf  in 

equation  (5)  is  63%  of  the  true  value  and  the  expected  value  of  A(t)  in 
equation  (3)  is  128%  of  the  true  value.  Trivially,  at  this  SNR,  the  expected 
value  of  the  received  power  is  twice  the  true  transmission  power  of  any 
slowly-varying  signal  (cf  equation  (4)). 

Signal  modulation  recognition  also  uses  time-averaged  sampled  characteristics 
such  as  the  variance  c2  of  the  signal  envelope  in  the  observation  period  T: 


o 


2 

A 


n 


A(j  t)  )  2 


where  t 


T/n. 


1  =  1 


(6) 


Consider  what  happens  to  this  variance  if  the  signal  strength  falls  so  that 
I B ( C ) |  is  significantly  less  than  o  for  the  entire  observation  period.  The 
approximation  E(A(t))  =  »'Va2/2  (1  +  B2(t)/4a2)  discussed  earlier  is  valid,  so, 
using  equation  (4) ,  the  expected  value  of  0^  is 


E(a2) 


£  \  (B2 (j t)  +  2a2)-  l  (2o2  +  |  )  (B2 ( j t )  +  (£  )  (B2(jt))2/8a2) 


(1  -  £)  (;j  ^  B2  ( j  t)  +  2c2) . 
j 


That  is,  the  expected  value  of  the  variance  of  the  envelope  is  proportional  to 
the  average  power  received.  It  is  independent  of  the  true  envelope  variance, 
and  so  carries  no  information  on  modulation  type. 

If  the  narrowband  noise  has  power  spectrum  symmetrically  distributed  about  the 

centre  frequency,  the  second  moment  of  0 *  ( t )  is  infinite.  (Although  this  is 

not  a  new  result,  it  is  not  well  documented  and  has  been  overlooked  in  at 

least  one  recent  paper(ref . 7) ,  so  a  proof  is  given  in  Appendix  I.)  However, 

any  least  squares  estimate  of  the  second  moment  is  of  course  finite.  If  the 

integrals  in  equation  Appendix  (1.7)  are  replaced  by  sampled  estimates,  the 

estimated  expected  value  of  (0,-2-rrf  )2  is  a  function  of  the  transmitted 

c 

signal's  envelope  B  and  its  time  derivative,  and  of  the  transmitted  signal's 
instantaneous  frequency  if'.  The  exact  dependence  is  complicated,  but  the 
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sensitivity  of  the  second  moment  to  transmitted  signal  parameters  is 
manifested  in  another  commonly-used  time-averaged  characteristic,  viz.,  the 
variance  of  the  signal  instantaneous  angular  frequency  over  the  observation 
period  T: 


o|,  =  (t/T)  \  0'2(jt)  -  ((t/T)  6'(jt))2 


To  illustrate  the  behaviour  of  this  characteristic  with  modulation  type  and 
noise,  various  signals  modulated  by  a  sinusoid  of  frequency  5  kHz  were 
simulated  with  the  procedure  described  in  Section  5.  The  same  Gaussian  white 
noise  was  added  to  each  signal,  the  total  power  was  normalised,  the  signal  was 
mixed  down  to  baseband  and  then  sampled  at  a  frequency  slightly  greater  than 
10  kHz  for  1024  samples.  Instantaneous  angular  frequency  was  calculated  as  in 
equation  (2).  This  was  repeated  at  various  signal  to  noise  ratios.  The  value 
of  Oq,  for  each  of  the  modulation  types  is  recorded  against  SNR  in  Table  1. 

The  theoretical  zero  variance  of  the  instantaneous  frequency  of  AM,  DSB,  and 
CW  signals  is  observed  only  at  100  dB.  Again,  the  difficulty  of  using  the 
variance  as  a  modulation  type  discriminator  without  reference  to  SNR  is 
apparent . 


TABLE  1.  STANDARD  DEVIATION  OF  0 '  BY  SNR  BY  MODULATION  TYPE  (x  10  kHz) 


SNR  (dB) 

AM 

DSB 

SSB 

CW 

FI 

*1 

(a) 

(b) 

(c) 

100 

0 

0 

0 

0 

■1 

13 

80 

0 

1 

0 

0 

1 

13 

40 

1 

127 

0 

0 

1 

13 

20 

9 

132 

3 

1 

2 

13 

10 

23 

72 

9 

4 

5 

14 

0 

29 

82 

42 

18 

18 

30 

—CO 

42 

42 

42 

42 

42 

42 

(a)  modulation  depth  0.8 


(b)  modulation  index  0.3 

(c)  modulation  index  3 


3.  SOME  ALTERNATIVE  SIGNAL  CHARACTERISTICS 

The  past  use  of  the  envelope,  frequency  and  zero  phase  as  basis  for  charac¬ 
terisation  of  signals  needs  no  explanation.  However,  it  is  clear  that  charac¬ 
teristics  which  provide  unbiased  estimates  in  noise  would  be  preferable,  as 
long  as  they  provide  as  good  discrimination  between  the  various  modulation 
types  in  the  absence  of  such  noise.  With  this  in  mind,  we  put  forward  some 
estimators  of  various  signal  characteristics  which  are  unbiased  in  the  pre¬ 
sence  of  narrowband  Gaussian  noise  symmetrically  distributed  about  the  centre 
frequency  f^. 
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In  the  following,  dependence  on  time  of  the  various  time  functions  is  not  made 
explicit.  An  overbar  on  a  time  function  denotes  the  sampled  time-average  over 
the  observation  period  T  ie: 

n 

X  =  |  V  X(jt),  T  =  nt  . 

j'l 

Define  cd  =  2frf  ,  8  =  <fc-U)  t  and  Y  =  0-to  t. 
c  c  c  c 

3 . 1  Lemma 

AA 

(i)  E(AA'  =  rr'+rr')  =  BB'. 

A  A  A 

(ii)  E(A2y'  =  rr'-r'r-co  (r2+r2))  =  B28'. 

c 

(iii)  If  the  noise  is  stationary  over  the  observation  period, 

E(A2-A2)  =  B2-B2. 

(iv)  If  the  noise  is  stationary  over  the  observation  period, 

E(A26'-A2e')  =  bV-bV. 

The  proof  of  (iii)  follows  trivially  from  equation  (A);  Civ)  follows  from 
parts  (ii)  and  (iii).  Appendix  II  outlines  the  proofs  of  (i)  and  (ii). 

Since  the  instantaneous  values  are  unbiased,  we  can  take  as  unbiased 

characteristics  of  the  signal  the  digitised  distribution  of  occurrences  of 
values  of  the  time  functions  A2y'  and  AA'  about  their  time-based  means. 

Note  that  we  are  here  talking  only  of  the  effects  of  noise,  not  of  the 
variation  between  signals  of  the  same  modulation  type  with  message. 

We  next  discuss  some  time  averaged  signal  parameters  which  could  provide 
useful  modulation  type  discriminators  and  which  are  unbiased  over  a 

sufficiently  large  number  of  observations,  as  long  as  the  noise  is 
stationary  and  the  noise  at  different  observations  is  uncorrelated  (a 
standard  assumption). 

Define 

n  n  n 

#(X, Y)  -  (2/n2)  7  X(1t)  )  Y (kt) -  i  V  X(jt)Y(jt)  =  2  X.Y  -  LY 

j=l  k-1  j-1 

That  is,  treating  X  and  Y  as  functions  of  time  #(X,Y)  is  the  product  of 
their  time-based  means  minus  their  time-based  covariance. 
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3.2  Lemma 

If  n  >>  1/SNR,  the  following  are  unbiased  estimators: 

(i)  #(a2,a2) 

(ii)  #(A2,AA') 

(ill)  //(A‘ ,  A2  (Y ’+c) ) ,  c  any  constant 
(iv)  //( AA' ,  A2  (Y'+c) )  ,  c  any  constant. 

The  proof  of  the  lemma  is  sketched  in  Appendix  II. 

Part  (i),  which  says  that  the  square  of  the  time  based  mean  of  the  squared 
envelope  minus  its  variance  is  an  unbiased  estimator,  does  not  need  the 

assumption  of  symmetry  in  the  noise  power  spectrum.  It  can  obviously  also 
be  expressed  in  terms  of  the  time-based  mean  and  variance  of  the  signal 
power. 

For  a  fixed  amplitude  signal  (FM,  CW  etc)  the  time-based  variance  of  B2  is 

zero,  so  that  the  expected  value  of  Section  3.2(i)  is  B2 ,  and  so  is 

proportional  to  the  square  of  the  transmitted  signal  power,  regardless  of 

the  total  noise-plus-signal  power  received.  This  parameter  is  potentially 

extremely  useful  even  with  amplitude-varying  signals.  For  instance,  if  a 

quartenary  AM  signal  is  observed  over  sufficient  interval  T  to  record  a 

more-or-less  equal  number  of  the  bit  levels  ±b  ,  ±2b  ,  then  the  time 

^  o  o 

average  of  B2  is  just  Kb2  +  (2b  )")  =  (5/2)b2  and  of  B4  is 

o  o  o 

i(b4  +  (2b  )4)  =  (17/2)b4.  Hence  #(BZ,B2)  =  4b4;  ie  the  expected  value  of 
o  o  o  o 

2  2 

the  observable  #( A  ,A  )  indicates  the  true  signal  amplitudes. 

4.  THE  SIGNAL  MODELS 

Consider  a  message  signal  m(t) ,  either  digital  or  voice,  transmitted  on  a 
carrier  of  frequency  f.  Write  the  transmitted  signal  S(t)  as 

S(t)  =  SI(t)cos(2^ft)  -  SQ(t)sin(2*ft).  (7) 

Table  2  lists  the  pair  (S^,Sq)  up  to  a  multiplicative  factor  according  to 
modulation  type:  the  lag  £  allows  for  asynchrony  between  message  and  carrier. 


7 
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TABLE  2.  QUADRATURE  PAIRS  BY  MODULATION  TYPE 


Modulation 

type 

SI 

SQ 

Notes 

AM 

l+am(t+Jl) 

0 

a  is  modulation  depth 

DSB 

m(  t+il) 

0 

SSB 

m(t+£) 

A. 

m(t+Jl) 

*  denotes  Hilbert  transform 

t+£ 

t+l 

FM 

cos(2iTb  m(x)dx) 

sin(2irb  /  m(x)dT) 

b  is  modulation  index 

‘0 

Jo 

FSK 

cos(m(t+ii,)A.  t) 

sin(m(t+Jl)  A.  t) 

A/2tt  is  frequency  jump 

PSK 

cos(rrm(t+jl)  /(a+1) ) 

sin(ir(m(t+£)  /  (a+1 ) ) 

data  set  {±1, . . .±a}. 

We  assume  that  the  transmitted  signal  essentially  lies  in  the  band  f  ±  B  : 

m 

that  is,  for  the  amplitude  modulated  signals  or  narrowband  FM  signals,  we 

require  m(t)  to  have  bandwidth  at  most  B  ,  whilst  for  wideband  FM,  we  require 

m 

this  bandwidth  to  be  much  less  than  the  modulation  index  b  where  b  £  B  .  For 

m 

the  digital  cases  PSK  and  FSK,  we  must  impose  conditions  on  the  pulse  shape  to 
meet  this  bandwidth  restriction.  With  S(t)  thus  constrained,  its  Hilbert 

A 

transform  S(t)  is  simply 


S(t)  =  Sq  cos(2irft)  +  Sj  sin(2iTft).  (8) 


Next,  suppose  the  incoming  signal  r(t)  is  received  in  the  band  f+d±B,  where  d 

is  the  offset  from  accurate  tuning  on  the  signal,  and  is  assumed  to  satisfy 

Id  I  <  B-B  so  that  the  transmitted  signal  lies  within  the  received  band.  The 
1  m 

incoming  signal  is  the  sum  of  the  transmitted  signal  S(t)  (passed  without 
distortion),  and  a  noise  component  n(t)  which  we  take  to  be  stationary, 
Gaussian  and  with  power  spectral  density  symmetrically  distributed  about  f+d. 
This  narrowband  noise  is  described  by 


n(t)  =  n[(t)  cos(2tr(f+d)t)  -  n2(t)  sin(2n(f+d)  t) 

*  Pi(t)  cos ( 2 7T f t )  -  p2(t)  sin(2irft) 


(9) 


where 

pt(t)  =  n[(t)  cos(2ndt)  +  n2(t)  sin(27rdt) 


and 


p2(t) 


-n2(t)  cos(2irdt)  +  nt(t)  sin(2irdt). 
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Its  Hilbert  transform  is  n(t)  =  P2(t)  .  cos(2irft)  +  pi(t)  .  sin(2mft).  The 
act  {n! (t) ,n2 ( t ) }  consists  of  observations  from  identically  distributed 
Gaussian  random  variables  with  variance  equal  to  the  average  power  of  S(t) 
divided  by  the  signal-to-noise  ratio  SNR.  With  a  rectangular  or  Gaussian 
bandpass  filter  acting  on  white  noise  we  can  suppose  the  observations  are  from 
independent  random  variables  when  sampled  at  the  Nyquist  rate. 

Finally,  we  suppose  that  in  obtaining  the  Hilbert  transforms  (cf  figure  1)  the 
signals  are  mixed  back  to  baseband  modulo  the  tuning  offset  and  a  constant 
phase  offset  6.  That  is,  we  set  f+d  *  0. 


5.  SIMULATING  THE  INFORMATION  SIGNALS  AND  NOISE 

Voice  message  signals  m(t)  were  obtained  via  a  64  kHz  digitisation  of  male 

t 

speech.  The  derivative  m'  and  the  integral  m  were  obtained  numerically  at 

‘  0 

this  sampling  rate.  The  resulting  data  were  subsequently  sampled  at  8  kHz  for 
AM  and  DSB  modulation,  and  32  kHz  for  wideband  FM.  An  example  of  this  voice 
data  is  shown  in  figure  2(a). 

Digital  message  signals  with  user  defined  symbol  duration  T  were  generated  by 

randomly  selecting  a  "data"  sequence  a0,  a; ,  a2  ...  from  the  set  {+ 1 , ±2 , . . . ±a] 

and  putting  this  information  on  a  Gaussian  pulse  of  bandwidth  B  truncated  to 

m 

a  3-symbol  duration.  Specifically,  setting  k  =  [t/T],  A  =  t-kT,  we  sampled 
m(t)  defined  for  kT  <  t  <(k+l)T,  k  >  1 ,  by 


m(t)  =  3(Vl 


-2ttB2  (A+T/2)  2 
m 


-2ttB2  (A-T/2)  2 

a  e  m  + 

k 


-2ttB2  (A+3T/2)  2 
m 


k+1 


).  (10) 


The  message  bandwidth  B^  and  sampling  rate  2B  were  user-selected;  to  avoid 

excessive  interference  between  adjacent  symbols,  only  B^T  products  greater 

than  l  were  selected.  Figure  2(b)  is  an  example  of  binary  message  data.  The 
sampled  derivative  m'  was  approximated  by  the  analytic  derivative  of 
equation  (10),  and  the  integral  of  m  was  obtained  numerically. 

The  Hilbert  transform  of  the  message  signals  was  not  computed,  so  that  SSB 
signals  were  not  simulated. 

The  message  signals  were  applied  to  the  various  modulation  types  as  in  §3,  to 
generate  the  transmitted  signals  S(t).  The  same  noise  n(t)  was  added  to  each, 
to  give  received  signals  r(t).  The  inphase  and  quadrature  noise  observations 
nj(t)  and  n2(t)  were  modelled  by  the  returns  from  a  Gaussian  random  number 
generator  with  variance  l/O+SNR)  where  SNR  is  the  user-supplied 
signal-to-noise  ratio.  The  average  received  signal  power  was  approximately 
set  to  1  by  normalising  S(t)  digitally  so  as  to  have  average  power 
SNR/(1+SNR).  For  example,  with  AM,  the  pair  )  of  Table  1  was  divided  by 
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n 


SNR  +  1 
2  n  SNR 


(i+a.m(t1))2)- , 


i=l 


where  m(t^)  was  the  sample  of  the  message  signal  ie  m(t^)  =  m((i+'.)/2B) 

for  i.  a  user-defined  lag.  The  number  of  samples  n  was  chosen  to  give  a  total 
observation  period  of  about  90  ms. 

The  Hilbert  transform  r(t)  was  readily  computed  from  the  equation 

r(t)  =  (SI(t)+p1  (t) )  cos(21Tdt+i))  -  (S  (t)+p2  (t) )  sin(2irdt+<£) .  (11) 


Since  the  functions  m' (t)  needed  to  compute 


S’(t) 


and  S^(t) 


were  supplied. 


the  derivatives  r'(t)  and  r'(t)  were  also  readily  computed,  given  the 
derivatives  nl(t)  and  nj(0  needed  to  compute  pl(t)  and  p|(t) .  These  were 
obtained  from  a  Gaussian  random  number  generator  whose  variance  depended  on 
the  bandpass  filter  modelled:  the  variance  was  tt2B2 / (6 fl+SNR) )  for  a 

rectangular  filter,  or  i~E2  /  ( 1+SN'R)  for  a  Gaussian  filter. 


6.  SIMULATION  RESULTS 

Figures  3  to  7  show  the  digitised  distributions  of  the  standard  variables  (the 
envelope  A  and  the  normalised  instantaneous  angular  frequency  9')  about  their 
mean  values,  plus  the  unbiased  distributions  of  A2,  AA'  and  A4?'  about  their 
means.  Each  plot  is  the  average  of  10  runs,  from  3  or  4  different  modulating 
signals  at  l  to  3  different  SNR  levels.  The  SNR  levels  were  in  the  range  15 
to  30  dE  for  strong  signals  (group  (a)  figures),  and  -1  to  3  dB  for  weak 
signals  (group  (b)  figures).  Group  A  figures  are  from  voice-modulated  AM,  DSB 
and  wide-  and  narrowband  FM  signals,  and  group  B  figures  are  from  binary  AM, 
DSB  and  narrowband  FM  signals.  The  corresponding  distribution  for  the 
Gaussian  white  noise  is  included  in  all  figures  for  comparison.  CW  signals 
are  Included  with  group  B  for  convenience  of  display.  The  marked  visual 
difference  between  group  (a)  and  (b)  figures  make  a  strong  heuristic  argument 
for  the  inclusion  into  any  classification  scheme  of  subclasses  related  to  the 
SNR. 

We  will  discuss  each  of  these  figures  in  turn,  but  caution  that  one  cannot 
fully  analyse  application  of  the  parameters  discussed  in  the  preceding 
sections  on  so  few  signals  and  modulating  forms.  Rather,  this  section 
indicates  the  relative  performance  of  the  various  parameters. 

In  the  following,  the  parameter  features  used  to  make  comparisons  are: 

(i)  the  number  of  occurrences  of  the  most  commonly  observed  quantised 
distance  from  the  mean:  that  is,  the  peak  of  the  distribution  and 
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(ii)  the  number  of  occurrences  of  distances  greater  or  less  than 
prespecified  limits:  that  is,  the  tails  of  the  distribution. 

The  2-sided  t-test.  is  used  to  determine  whether  mean  values  recorded  for 
the  various  modulation  types  are  the  same  within  a  99%  confidence  level. 
(Here,  the  mean  value  referred  to  is  the  average  value  of  the  histogram  bin 
recorded  over  the  10  simulation  runs.  In  using  this  test,  we  are 
implicitly  assuming  that  the  number  of  observations  which  lie  in  a 
particular  histogram  bin  is  a  normally  distributed  random  variable  with 
variance  independent  of  the  modulation  type.  However,  the  t-test  is  known 
to  be  robust.)  We  use  standard,  albeit  statistically  loose,  terminology 
below  in  saying  that  "two  modulation  types  can  be  distinguished"  when  the 
t-test  fails,  and  by  saying  "the  values  are  the  same"  if  it  is  passed. 

(1)  Figure  3:  Distribution  of  envelope  A  about  mean 

In  figure  3,  the  envelope  of  the  noise  approximates  the  Rayleigh 
distribution  which  in  the  large  sample  limit  it  must  follow.  The  envelope 
of  a  DSB  signal  is  the  absolute  value  of  the  modulating  signal.  As  a 
consequence  of  the  randomness  essential  to  any  information-carrying  signal, 
the  distribution  of  the  voice  DSB  envelope  is  very  similar  to  the  Rayleigh 
distribution  for  the  voice  signal  (figure  3A(a)).  The  envelope  of  the  AM 
signal  bears  the  influence  not  only  of  the  message  but  also  of  the  carrier, 
so  for  strong  signals  has  a  higher  peak  at  the  mean.  Finally,  the  FM  and 
CW  strong  transmissions,  with  their  almost  constant  envelopes,  have  sharply 
spiked  distributions. 

However,  for  weak  signals,  ALL  the  modulation  types  have  the  same 
distribution  of  envelope  as  has  the  white  noise.  (Figure  3A(b)  and  3B(b)). 
Modulation  discrimination  is  unreliable  using  these  distributions,  as 
foreshadowed  in  Section  3. 

(2)  Figure  4:  Distribution  of  instantaneous  angular  frequency  8'  about 
mean 

In  figure  4,  the  distribution  for  noise  of  the  instantaneous  frequency 
0 ’ / 2tt  about  its  mean  approximates  a  density  function  of  the  form 

a(l+0,2/b)  (see  Appendix  III).  For  strong  voice  signals 
(figure  4A(a)),  the  narrowband  FM  distribution  follows  a  noise-like 
distribution,  again  reflecting  the  randomness  of  the  modulating  signal. 

In  line  with  the  sort  of  results  recorded  in  Table  1  and  discussed  in 
Section  3,  the  distribution  for  strong  CW  signals  appears  sharply  peaked, 
but  approaches  that  of  the  noise  when  the  signal  is  weak  (figure  4B) . 
Since  the  strong  signal  category  includes  SNR's  down  to  15  dB,  there  is 
sufficient  noise  to  broadly  spread  the  mean  DSB  distributions 
(figures  4A(a)  and  4B(a)),  although  as  with  the  CW  signals  there  is  a  large 
variance  associated  to  each  bin  value.  On  the  10  runs  performed  on  strong 
CW  signals,  the  mean  histogram  values  returned  were  not  significantly 
different  from  zero  -  with  a  one-sided  t-test,  99%  confidence  level  - 
or  from  the  corresponding  strong  voice  FM  histogram  values.  Although 
figure  4A(a)  suggests  a  moderately  peaked  distribution  for  the  strong  AM 
voice  signals,  none  of  the  histogram  values  were  in  fact  signif icantly 
different  from  the  corresponding  white  noise,  weak  FM,  CW  or  strong  digital 
FM  values.  The  remaining  AM  and  DSB  signal  types  (strong  or  weak,  voice  or 
binary)  could  not  be  distinguished.  Only  the  strong  voice,  FM  signals 
could  be  separated  from  the  various  FM  signal  types. 
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Figure  3.  Distribution  of  envelope  A  about  mean 


Figure  4.  Distribution  of  instantaneous 


Figure  5.  Distribution  of  Az  about  mean 
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(3)  Figure  5:  Distribution  of  A2  about  mean 

In  the  large  sample  limit,  the  distribution  of  the  square  of  the  envelope 
of  white  noise  is  the  exponential  distribution(ref . 5) .  Shifted  by  the 
mean,  it  has  asymptote  at  -2,  because  the  average  value  of  the  envelope 
squared  is  twice  the  average  receiver  power  which  we  have  normalised  to 
unity.  Comments  made  in  Subsection  (1)  concerning  the  envelope  generally 
apply  here,  with  the  important  exception  that  for  noisy  signals,  only  the 
digital  DSB  cannot  be  statistically  distinguished  from  noise. 

(4)  Figure  6:  Distribution  of  AA'  about  mean 

The  variable  AA'  is  proportional  to  the  derivative  of  the  envelope  squared. 
For  the  white  noise  considered,  it  is  the  sum  of  products  of  independent 
Gaussian  random  processes  (cf.  3 . 1 ( i ) )  and  so  is  itself  normally 
distributed.  This  is  observed  in  figure  6. 

Visually,  three  groups  can  be  recognised  for  the  strong  signals  -  AM/DSB, 
wideband  FM/noise,  and  narrowband  FM/CW  (figures  6A(a)  and  6B(a)). 
However,  statistically  the  FM,  CW  and  noise  distributions  coincide  and 
voice  AM,  binary  AM,  and  DSB  are  mutually  distinguishable.  The  addition  of 
noise  (figure  6(b))  does  not  significantly  alter  the  distributions  of 
digital  DSB  and  AM  but  for  the  voice  data  spreads  the  AM  distribution  to 
make  it  significantly  different  to  all  other  signal  types  considered.  This 
is  also  the  case  for  the  wideband  FM  distribution.  The  remaining  FM  signal 
types,  CW  and  noise  could  not  be  statistically  distinguished. 

(5)  Figure  7:  Distribution  of  A26'  about  mean 

Following  equations  (2),  the  variable  A26'  can  be  expressed  as  a  difference 
of  products  which,  for  the  noise  signal  considered  here,  are  Gaussian 
random  processes.  Hence,  as  observed  in  figure  7,  the  noise  distribution 
is  approximately  normal.  For  the  strong  CW  signal,  it  is  sharply  peaked 
(figure  7B(a) ) . 

Comments  made  about  the  distribution  of  instantaneous  frequency 
(Subsection  (2))  generally  apply,  except,  predictably,  more  discrimination 
between  the  various  AM  and  DSB  signal  types  is  possible. 

As  a  comparative  summary  of  the  results  obtained,  Tables  3  and  4  display 
symmetric  binary  matrices,  where  a  zero  in  the  (i,j)th  position  indicates  that 
all  features  (mean  histogram  values)  in  the  set  under  consideration  were  the 
same  for  the  ith  and  jth  signal  types  (t-test,  99%  confidence  limit).  The 
features  are,  for  Table  3,  the  peak  and  tails  values  of  the  histograms  of  the 
envelope  A  and  instantaneous  angular  frequency  6'  and,  for  Table  4,  the  peak 
and  left  tail  value  of  the  histograms  of  A2,  A.',  and  A29'  (a  total  of  six 
features  in  each  case).  Using  the  standard  parameters  there  are  14  unresolved 
signal  pairs,  whilst  with  the  new  parameters  there  are  only  two. 

Table  5  lists,  for  the  various  modulation  types  and  for  strong  and  weak 
signals,  the  mean  and  the  associated  standard  deviation  over  the  10  simulation 
runs  of  the  variable  #( A2, A2)  defined  in  Section  3.2(i).  As  expected,  this 
parameter  appears  to  be  a  very  good  way  of  determining  the  true  signal  power, 
and  hence  the  SNR,  of  constant  envelope  signals.  The  other  variables  defined 
in  Section  3.2  proved  to  have  large  standard  deviations  within  each  signal 
type,  and  so  appear  to  have  limited  use  in  modulation  recognition.  However, 
this  very  variability  may  make  them  useful  fine  tuners,  for  instance,  in 
distinguishing  between  voice  and  digital  modulating  signals,  given  an  estimate 
of  SNR. 


ERL-0367-TR 


18 


TABLE  3.  SIGNIFICANTLY  DIFFERENT  FEATURE  VALUES  BY  SIGNAL  TYPE  (A, 6') 
(t-TEST,  99 %  CONFIDENCE  LEVEL) 

Features:  peak  and  tails  of  distributions  of  envelope  A  and 
instantaneous  angular  frequency  9'. 
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Strong  binary-Ba 
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AM 

DSB 

FM 

FM 

AM 
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AM 
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CW 

AM 
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CW 

Noise 
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nb 
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nb 

nb 
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0 

1 

i 

0 

1 

1 

1 

1 

1 

1 

i 

1 

fl 

i 

1 

1 
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1 
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TABLE  4.  SIGNIFICANTLY  DIFFERENT  FEATURE  VALUES  BY  SIGNAL  TYPE 
(A2,  AA\  A26 ' )  (t-TEST,  99%  CONFIDENCE  LEVEL) 

Features:  peak  and  left  tail  of  distributions  of  A2, 
AA'  and  A2e' 
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nb  -  narrowband 
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TABLE  5.  THE  MEAN-SQUARED-ENVELOPE  SQUARED  MINUS  ITS  VARIANCE  If  ( A2, A2) 


AM 

DSB 

Narrowband 

FM 

Wideband 

FM 

CW 

Noise 

(a) 

(b) 

(a) 

(b) 

(a) 

(b) 

(a) 

Strong  signals 

■ 

mean 

3.1 

1.0 

-4.0 

BBS 

3.9 

3.9 

3.9 

3.8 

sd 

0.3 

0.7 

1.4 

HI 

0.1 

0.1 

0.1 

0.1 

EH 

Weak  signals 

■ 

■ 
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n 

-1.1 

-0.4 

H 

1.1 

1.2 

1.1 

sd 

i 

EH 

0.7 

0.6 

0.4 

0.4 

0.4 

_ 

(a)  voice  modulated  signal; 

(b)  digitally  modulated  signal 


7.  CONCLUSIONS 

We  have  shown  theoretically  why  features  based  on  the  standard  time  domain 
parameters  of  signal  envelope  and  instantaneous  frequency  are  strongly 
influenced  not  only  by  the  modulation  type  but  also  by  the  signal  to  noise 
ratio.  Parameters  which  are  unbiased  in  the  presence  of  Gaussian  white  noise 
have  been  introduced  and  limited  simulations  run  to  compare  their  performances 
with  those  of  the  standard  parameters. 

Discrimination  between  transmissions  of  modulation  type  AM,  FM,  DSB  or  CW  on 
the  basis  of  characteristics  of  the  distribution  of  the  new  parameters  appears 
to  be  at  least  as  good  on  strong  signals,  while  in  a  noisy  environment 
modulation  recognition  is  markedly  enhanced.  The  further  discrimination 
between  voice  and  binary  data  using  these  techniques  may  also  be  feasible. 
More  work  on  the  statistics  of  the  speech  and  multilevel  digital  modulating 
signals  needs  to  be  done.  Performance  in  non-Gaussian  noise  also  needs  to  be 
evaluated  both  theoretically  and  practically. 

Further  simulations,  on  a  much  larger  data  base,  and  extensions  to  SSB,  00K, 
PSK  and  FSK  are  planned  or  under  way.  A  hardware  preprocessor  is  being  built 
to  provide  analytic  computation  of  the  new  variables,  which,  after 
digitisation,  will  aid  in  providing  a  sufficiently  large  signal  set  for  the 
implementation  of  a  classification  scheme.  It  is  envisaged  that  the  Fisher 
linear  discriminant  function  will  be  applied  to  features  obtained  from  the 
parameter  distributions. 
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APPENDIX  I 

VARIANCE  OF  THE  INSTANTANEOUS  FREQUENCY 

Suppose  a  received  signal  r(t)  «  A(t)cos(0(t) )  is  the  sum  of  a  transmitted 

signal  S(t)  -  B(t)cos($(t))  and  narrowband  Gaussian  noise  n(t)  =  ni(t)cos(o)  t) 

c 

-  n2 (t) sin(a)ct)  which  has  power  spectrum  symmetrically  distributed  about  the 
centre  frequency  w  /2ir.  Then  the  expected  value  of  6' 2  is  not  defined  (is 
infinite) . 


Proof.  Let  y  “  0-ui  t  and  6  •  dr-ui  t. 
-  c  c 

has  inphase  component  R^  •  BcosB+nj, 
It  is  straightforward  to  compute 


With  respect  to  the  centre  frequency,  r 
and  quadrature  component  R  =  Bsin8+n2. 


r'  =  R'  cosuj  t  -  R’  sinto  t  -  to  r 
I  c  Q  c  c 

and 


r'  =  R'  coso)  t  +  R'  sinw  t  +  a)  r. 
Q  c  I  c  c 


Moreover  (cf  equation  (2)), 


y'  -  0 ' -w  =  (r(r'-o)  r)  -  r(r'+co  r))/(r2+r2). 

c  c  c 


From  equations  (1.1)  and  (1.2),  we  get 


(^-RIRQ)/(R!+RQ>- 


(1.1) 


(1.2) 


(1.3) 


Substituting  the  definitions  of 


Rj  and 


Rp  gives 


y'  “  { (Bcos6+ni)  (B'sin8+6'Bcos8+n2)  -  (B '  cosg-B '  Bsin8+n  } )  (Bsin£?+n2 )  } 

/{ (BcosB+n! ) 2  +  (Bsin£+n2) 2 }.  (1.4) 


v 


Recall  that,  at  any  instant  t,  n  i,  n  2,  n'i  and  n'2  are  independent  random  vari¬ 
ables  with  zero  mean.  Moreover,  E(nf)  -  E(n2)  -  o2,  and  E(n!2)  *  E(n22  ). 
(This  is  a  consequence  of  our  assumptions  on  the  form  of  noise.)  It  is  then 
easy  to  show  from  equation  (1.4)  that 
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E ( y * 2 )  =  E({(B'sinB+6'Bcos8) (Bcosg+nj)  -  (B* cosg-B' Bsing) (BsinB+n2) }2 

/{  (Bcosg+ni )  2  +  (Bsing+n2)2}2) 

+  E(nJ  2) .  E(l/{  (Bcosg+ii! )  2  +  (Bsing+n2)  2  }) .  (1.5) 


Define  a  =  y~B>  so  Chat  nj+Bcosg  =  Acos(a+8)  and  n2+Bsing  =  Asin(a+g).  A 
short  calculation  then  reduces  equation  (1.5)  to 


E(y'2)  =  E((-B'sina+6'Bcosa)2/A2)  +  E(nl2)  E(l/A2).  (1.6) 


Now  because  n  1  and  n2  are  independent  Gaussian  random  variables,  they  have 
joint  density  function 

p(n[n2)  *=  (1/2t to2)  1  exp (-(nf+n|) /2o2) )  =  (1/2-jto2)  1  exp(-(A2+B2-2ABcosa) /2o2) . 

(1.7) 


From  the  Jacobian  calculation,  we  have  dn:dn2  =  AdAda.  Hence  equation  (1.6) 
becomes 


E(y' 2 )  =  (1/2’tc:)  1  e 


-B2/2 a2 


da  ((-B'sina+B'Bcosa)2  +  E(n'i2)) 


CO 

dA(e"(A2'2ABcOS&)/2o2))A).  (I’8) 

0 


Set  1.(6)  =  dA(e  ^A  +^AB)/2a  Clearly  1(6)  does  not  converge  as 

A  A 

'  6 

as  6-*0.  But  from  equation  (1.8),  if  620, 

.2^ 

E  ( y ' 2 )  >  (e"B2/2°  /  (2rc2) )  I.  (6)  (  da(-B '  sincd-g'  Bcosa) 2  +  2tt  E(n',2)). 

A  / 

'  o 


This  implies  our  result. 
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APPENDIX  II 

SOME  UNBIASED  ESTIMATORS 

To  prove  that  the  functions  defined  in  Lemmas  3.1  and  3.2  are  unbiased 
estimators  as  claimed,  one  can  proceed  as  follows. 

We  have 


r(t)  =  S(t)  +  n(t) , 

r (t)  =  S(t)  +  n(t)  (II. 1) 


whe 


re  r  =  Acos0  =  Acos  (y+^t)  and  S  =  Bcos^j  =  Bcos  (B+Wct)  • 


The  noise  n  has  input  power  o2*  If  n^.  and  n^  are  respectively  its  inphase  and 
quadrature  components  with  respect  to  the  centre  frequency  f^,  then  n^.n^  and 
their  time  derivatives  n^  and  n^  are  (at  any  instant  t)  independent  zero-mean 
Gaussian  random  variables.  The  variance  of  n^  and  n^  is  a2  •  and  the  variance 
of  n^  and  n^  is  C(B)02,  where  the  form  of  the  function  C  of  noise  bandwidth  B 
is  determined  by  the  noise  power  spectrum  shape.  Now 


n'  =  nlcosat  t  -  n'sin^  t  -  w  n 
I  c  Q  c  c 

n'  =  n'cosco  t  +  n'sinoj  t  +  u  n.  (II. 2) 

y  c  i  c  c 


Using  equation  (II.  2),  and  the  independence  of  the  random  variables  n^,  n|,  n^ 


and  n',  it  is  straightforward  to  compute 


E  (nn) 

=  E  (n  (cos  2oj( 

A  A 

E(nn’ ) 

=  0  =  E(nn' ) 

E(nn' ) 

=  -u,cE(n2)  = 

E(nn'  ) 

=  wco2; 

E  (  n '  2  ) 

=  C(B)o2  +  ix,co 

I  Q' 


=  0; 


"Wco‘ 


(II. 3) 


Finally,  we  will  need  the  fact  that  if  Uj  ,  u2  ,  u3  and  u4  are  zero  mean 
,  Gaussian  random  variables,  then  E(u;U2u3)  =  0  and  E(ujU2u3u4)  *  E(uiU2)E(u3u4) 

+  E(u,u3)E(u2u4)  +  E(u jU4) E(u2u3) ,  (reference  5,  page  97). 
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To  prove  Section  3.1(11)  viz;  E(A2y')  -  B2  B  1  ,  use  the  definitions  in 
equation  (2) ,  to  write 

A  A  A  A  A  A 

A2v'  =  (rr'-rr')  -  oi  (r2+r2)  =  r(r'-oj  r)  -  r(r'-Hjj  r) . 

c  c  c 

Rewriting  as  in  equation  (II.  1),  this  is 

A.  A.  A  A  A  A 

(S+n)  (S'+n'-u)  (S+n))  -  (S+n)  (S '+n'+u)c (S+n) ) 

=  S ( S '  —<jJc S )  -  S(S'-hccS)  +  nn'  -  u^n2  -  nn'  -  a^n2  +  terms  linear  in  noise. 


Using  equation  (II. 3)  E(nn'-nn')  =  cjcE(n2+n2),  and  of  course,  the  expected 
value  of  any  term  linear  in  noise  (ie  in  n,  n,  n'  or  n')  is  zero. 

Thus 


E(A2y')  =  S(S'-ucS)  -  S(S'-cocS)  =  B2B’  . 

This  proves  Lemma  3. 1 (ii) -  Part  (i)  of  Lemma  3.1  is  proved  similarly. 

The  analysis  of  the  time-averaged  quartic  expressions  of  Section  3.2  is  more 
tedious,  but  just  as  elementary.  As  instance,  we  will  prove  part  (iii)  with 

c=0 ,  viz.,  E ( 2  A2.AV  -  AV)  =  2  B2.B2B’  -  B4B'. 

Using  equations  (2)  and  the  definition  y  =  6-w  t, 

A  A  A  A 

A4y'  =  (r2+r2)  (r(r'-0)cr)  -  rtr'-k^r))  =  fi(r)-f2(r), 

where  for  notational  convenience  we  define 

f  i  (r)  =  (rz+  r 2 )  (r (r ' -co^r ) ) ;  f2(r)  =  (r2+  r2) (r (r '+  wcr) ) . 
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f  l  (r)  »  (S  +2Sn->-n  +S  +2Sn+n  )  (SS'+Sn’+nS'+nn'-co  (S  +2Sn+n  )) 

c 

A  A  A  A  A  A 

=  (S2+S2)  (S(S'-U)  S))  +  (S2+S2) (nn'-w  n2)  +  (n2+n2)  S(S'-w  S) 

c  c  c 

AA  A  A  A  A 

+  (2Sn+2Sn)  (n(S'-20JcG)+n'S)  +  (n2+n2)  n(n'-U)  n) 

+  terms  linear  or  cubic  in  noise. 


Using  equations  (II. 2)  and  (2),  we  get 


E(f l (r) ) 


(S  +S  )  (S(S'-O)  S)  +  (<d  cr-io  a2))  +  2o2  S(S'-u)  S) 
c  c  c  c 


+  2 S ( S  ’  —  2uj^ S ) ct 2  +  2S2Wc<j2  +  4g2  (u^o2-u}^o2) 

A  A  A  A 

=  (S2+S2)  S(S’-u)  S)  +  4o2  S(S'-O)  S)  =  f  1  (S)  +  4o2  S(S'-oj  S), 

c  c  c 


Similarly, 


E(fz(r))  =  f 2  (S)  +  4 o2  S(S'-KJ  S) . 

c 


E(AV)  =  E(fi(r)-f2(r)) 


(II. 4) 


fi(S)  -  f z  (S)  +  4G2 (S(S'-w  S)  -  S (S '-ho  S))  =  B’6'  +  4a2B26'. 

c  c 


Recall  that,  given  n  samples  X(t) , . . . .X(nt)  of  a  random  process  X,  we  have 


defined  X  =  (1/n)  )  X(jt).  If  X  and  Y  are  uncorrelated  at  different  sample 


times. 
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E(X.Y)  =  (1/n2)  E(X(jt))E(Y(kt))  +  ^  E(X(  j  t)  Y(  j  t)  )  ) 

j 

=  EOO.E(yT  +  i  (eT3cyT  -  E(X)E(Y) )  .  (II. 5) 

Equation  (II. 5)  is  valid  for  the  case  X  *  A2 ,  Y  =  A2y'  because  of  our 

assumptions  about  the  noise.  With  equations  (4),  (II. 4)  and  (II. 5),  and 

Lemma  3.1(i),  we  have 

E(AV-2A2.AV)  =  E(AV)  “  2E(A2)  E(A2y»)  -  (2/n)  (EfA-y')  -  E(A2 ) E(A2y ' ) ) 

=  B46'+  4o2B26'  -  2 (B2+2o2 )  B2e'  -  (2/n)  (2C2  B26') 

=  B-g'  -  2B2 . B2 8 '  -  (4/n)  (o2B26'). 

We  can  neglect  the  last  term  if  2B2  >>  (4/n)a2,  ie  if  n  >>  1/SNR. 
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APPENDIX  III 

PROBABILITY  DISTRIBUTION  OF  THE  INSTANTANEOUS  FREQUENCY 

Assume  a  Gaussian  white  noise  signal  with  a  band-limited  power  spectrum 
symmetrically  distributed  about  its  centre  frequency  and  quadrature  components 
r.j  =  rcos6  and  n^  «  rsin0.  Let  p(0*)  be  the  probability  density  function  of 

the  instantaneous  angular  frequency  of  this  signal  about  its  centre  frequency. 

Suppose  o2  is  the  power  of  the  noise,  and  t*2  that  of  its  time  derivative. 

Then  (reference  7), 


p(S')  =  p(r,0,r',0')  dbdrdr ' , 


where  p(r , 0 , r ' , G ' )  can  be  determined  from 


P(nQ’ni,nQ,ni) 


)  exp(-(n2+n‘ ) /2o2 )  exp (- (n^2+n j 2 ) /2e2 ) 


Specifically  (cf  equation  (1.7)  to  (1.8)), 


p('r')  =  (4  tc2e2)  1  r2  exp  (-r  “ /2c2 ) .  exp  ( (r '  2+r  20  1 2 ) /2<?2 ) .  d@ .  dr .  dr ' , 

•  CJ0 


so  that 


p(e  ’ )  «  r2  exp(-r2(  1+0' 2 .  '2  /**2 )/ 2c2 )  .dr 


( 1+0' 2 . (a2 /«2 ) ) 


Note  that  c2 /&2  is  a  function  of  the  power  spectrum  shape  and  varies  as  B  2 
with  the  bandwidth  B. 
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